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The impact of quantum fluctuations on the phase diffusion in resistively shunted superconducting 
quantum points subject to an external ac-voltage is studied. Based on an extension of the classical 
Smoluchowski equation to the quantum regime, numerical results for fractional and integer Shapiro 
resonances are investigated to reveal characteristic features of quantum effects. It is shown that 
typically in the quantum regime the broadening of resonances cannot be described simply by an 
effective temperature. 



I. INTRODUCTION 

In recent years there has been a notable advance 
in the understanding of electronic transport through 
superconducting nanosystems. In particular, the de- 
velopment of fabrication techniques such as scanning 
tunneling microscopy, break-junction and lithographic 
methodologies 1 - have allowed to study atomic-size metal- 
lic contacts as ideal systems to test fundamental prop- 
erties of charge transfer through superconducting weak 
links. These achievements have not only deepened our 
understanding of subgap structures in the current-voltage 
characteristics but also revealed even microscopic de- 
tails of the contact such as transmission coefficients 2 ^. 
Theoretical predictions based on Greens-function or 
mean-field approaches have been confirmed in various 
experiments^. In this context, set-ups where atomic- 
size tunnel junctions are subject to both dc- and ac- 
voltages give access to the intimate relation between 
phase dynamics, driving, and dissipation leading to pro- 
nounced Shapiro resonances of integer and fractional or- 
der. 

Similar to conventional weak links, atomic point con- 
tacts are characterized by two energy scales^, namely, 
the coupling energy between the superconducting do- 
mains (Joscphson energy) and the charging energy of 
the junction. The competition between these two scales 
is crucially influenced by the electromagnetic environ- 
ment so that a realistic modeling of charge transport 
across the contact must necessarily incorporate its em- 
bedding in an actual circuit. Accordingly, the dynam- 
ics of the superconducting phase difference as the only 
relevant degree of freedom exhibits a diffusive motion 
subject to noise which, as it is well-known from the 
physics of Josephson junctions*, can be visualized as 
the Brownian motion of a fictitious particle. It turned 
out that for atomic point contacts this phase dynam- 
ics occurs in an overdamped regime. The corresponding 
classical frame is provided by the Smoluchowski equa- 
tion which has already been the starting point for cal- 



culations of current-voltage characteristics of Josephson 
junctions in low impedance environment s 10 ' 1 1 . There, 
only the Joscphson energy remains as relevant parameter 
while charging effects related to inertia drop out. The im- 
pact of quantum fluctuations have been attacked within a 
time-dependent perturbation theory ir>i£ where the whole 
range from coherent to incoherent Cooper pair transfer 
in the domain of Coulomb blockade could be captured. 
Later, a generalization of the Smoluchowski approach to 
the quantum regime (Quantum Smoluchowski) developed 
by one of us (JA) and co-workers^— has allowed to de- 
rive the same physics in a very elegant manner and in 
close analogy to the classical description. Quantum noise 
has been shown to be inevitably associated with charging 
effects according to the uncertainty principle, thus phys- 
ically ruling the changeover toward Coulomb blockade 
dominated transport. 

The motivation for the present work is two-fold. On 
the one hand, it is based on experiments conducted with 
atomic point contacts in the last years in the Quantron- 
ics group 12 ' 14 and on the other hand on a correspond- 
ing description in terms of the classical Smoluchowski 
approach^. While experimental results with ac-driven 
junctions followed theoretical predictions for the struc- 
ture of Shapiro resonances, substantial discrepancies ap- 
peared for their heights and widths. A plausible expla- 
nation has been the presence of residual spurious noise 
sources which lead to an effective temperature at the con- 
tact different from the actual base temperature. Here, we 
analyze if and if yes to what extent quantum noise must 
also be incorporated into this picture. Eventually, this 
may open the door to unambiguously characterize quan- 
tum effects in overdamped systems at low temperatures. 

The paper is organized as follows. In Sec. [IT] we present 
our generalization of the RS J model for contacts of arbi- 
trary transmission in the presence of microwave radiation 
and in the presence of quantum fluctuations. In Sec. IIIII 
the numerical method to solve the generalized quantum 
Smoluchowski equation is outlined together with the rel- 
evant scales and approximations. Section lTVl discusses re- 
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suits for the I —V characteristics and fractional Shapiro 
resonances, before in Sec. [V] the question whether quan- 
tum noise can be captured by an effective temperature 
is addressed. Section IVII is devoted to discussion and 
conclusions on future experimental realizations. 



II. THE MODEL 

We model a superconducting tunnel junction using 
an equivalent circuit with so-called lumped circuit pa- 
rameters that includes both the effect of dissipative 
sources and the distributed capacity. For weak links 
the standard resistivcly and capacitivcly shunted junc- 
tion (RCSJ) model captures the essential physics even 
in presence of an external ac-voltage^. The equivalent 
circuit, shown in Fig. [TJ is formed by a contact with a 
resistance i?, a capacitance C, and biased by a dc-current 
lb- The superconducting phase difference across the junc- 
tion denoted by 9 is the only relevant degree of freedom 
with a current-phase relation 1(0). Energy dissipation 
in the contact is accompanied by Johnson-Nyquist noise 
I n (t) thus the current conservation relation gets: 
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FIG. 1: Lumped circuit model of a superconducting weak 
link with capacitance C and within a resistive environment R 
subject to a dc-bias current h and an ac-voltage. 
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Here, the first term on the RHS is the displacement 
and the third term is the dissipative current. Further, 
V to t(t) = V(t) + V ac cos(ujt) with voltage V across the 
contact and V ac being an additional ac-voltage induced 
by a microwave field. According to V — $0$ ($o = fi/2e) 
this equation is in fact an equation of motion for the 
phase, i.e., 

^2 + I b = $oC0(t) + !£d(t) + 1(9) + I n (t) , (2) 

which is equivalent to the Brownian motion of a ficti- 
tious particle with mass M = &qC, friction constant 



r\ = 1/RC and classical noise force Zci(t) = —RI n (t)/$o 
in the potential 



U(6,t) = $ / I(9)d9- $ Q 9[I b 
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The noise force has zero mean and obeys 
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with D C i = k B T= 1//3. 

The regime where the capacitance is negligible thus 
corresponds to the strong friction domain (Smoluchowski 
regime) in the mechanical analog^. Previous work has 
shown that this is indeed the range where phase diffu- 
sion in superconducting atomic point contacts happens 
to occu r 12 ' 13 . Strong friction considerably simplifies the 
description and for the situation with zero driving even 
allows for analytical results. It is then often convenient to 
switch from the classical Langevin equation correspond- 
ing to ©, i.e., 



1 dU(9,t) 
r]M d9 



+ Z C i(t) 



(5) 



to an equation of motion for the probability distribution 
P(0,t), namely, 



dP(9,t) 
dt 



1 d 
^Md9 



dU(0,t) d 
'^9— + Dci d9 



P(0,t). (6) 



As first pointed out ir>i£, this classical Smoluchowski 
equation (SE) can be generalized to the low tem- 
perature domain where quantum fluctuations become 
substantial^. The quantum Smoluchowski equation 
(QSE) has been studied since then in a variety of 
applications 2 ^, including particularly an extension of the 
classical Ivanchenko-Zil'berman theory for Josephson 
junctions in low impedance environment o 10 ' 11 . There, 
quantum fluctuations are related to charging effects and 
reveal signatures of Coulomb blockade physics. The 
QSE follows from its classical counterpart by replacing 
Dei — > Dq(9) with the position dependent quantum 
diffusion coefficient 



Dq(9) 



1 - A/3U"(9) 

with a friction and temperature dependent function 



A = 2p 



2n 2 p 
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where \& denotes the digamma function and c = 
0.5772 . . . Euler's constant. Further, the charging en- 
ergy is E c = 2e 2 /C and we introduced the dimension- 
less resistance p = R/Rq with Rq = h/Ae 2 . Usu- 
ally, p <C 1 for circuits operated in the overdamped 
regime. The classical Smoluchowski range corresponds 
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to the high temperature limit r)h(3 = (3E c /(np) 1, 
where A w f3E c /ir 2 <C 1, while at low temperatures 
pE c /(jrp) 3> 1 quantum fluctuations are substantial ac- 
cording to A rj 2p\n(f3E c /ir 2 p). The generalization of 
the classical Langevin equation follows from ([5]) by re- 
placing Zci Zq = y/P Dq(8) Zci which describes a 
classical stochastic process with multiplicative noise. 

In the sequel we consider an atomic point contact 
with one conduction channel with transmission proba- 
bility t € [0,1]. Generalizations are straightforward. 
As it is well-known the current through the contact is 
then carried by two Andreev bound states with energies 



E±(0,t) = ±Av/l - Tsin 2 (6>/2) (A is the superconduct- 
ing gap) . If we restrict ourselves to voltages much smaller 
than A and fceT, there are no Landau-Zener transitions 
between Andreev states and an adiabatic approximation 
applies. Thus, the current-phase relation gets 
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This expression simplifies to the known sinusoidal rela- 
tion for tunnel junctions in the low transmission limit 
(t — > 0) and is proportional to sin(0/2) in the ballistic 
limit for r — > 1. In this latter domain and for exter- 
nally driven contacts higher harmonics become relevant 
such that apart from the conventional integer Shapiro 
steps also fractional ones appear. This situation has been 
studied in the classical realm in RefJ^. Here, we focus on 
the low temperature region where the classical descrip- 
tion must be extended to include quantum fluctuations 
as discussed above. 

Before we proceed, let us specify the domain in pa- 
rameter space where the strong friction approach and the 
modeling of the environment in terms of an ohmic resistor 
apply. With respect to the first issue we consider circuits 
of the type shown in Fig. [TJ Then, roughly speaking, the 
friction constant must sufficiently exceed all other rele- 
vant frequencies. In the mechanical analog this means 
r\ ^> ujjh/3, oJj/rj, eV ac /H, w with plasma frequency ujj 



\fEjE c /h and Josephson energy Ej = A(l — y/l — r). 
Note that the last condition rj 3> u> enures that the ex- 
ternal driving acts on time scales sufficiently larger than 
the relaxation time for momentum which is of order 1/r/. 
In terms of circuit parameters one has 

> npEj,(3E c E j7 eV ac , thj (10) 

irp 

with p <C 1. In addition, as discussed above, the ratio 
Hf3r) = f3E c /(irp) controls the impact of quantum fluctu- 
ations. Further, the adiabatic description ^ is justified 
if Huj < 2AV1 - t to avoid driving induced mixing of 
the two Andreev surfaces. 

With respect to the second issue, the modeling of the 
environment as being purely ohmic is of course a crude 
approximation to actual experimental set-ups. Any real- 
istic circuit exhibits at least a cut-off frequency fl c due to 



unavoidable additional capacitances. The Smoluchowski 
description remains valid as long as there is still a time 
scale separation between relaxation in phase [approach 
of a quasi-stationary state for P(0, t)] and the response 
time of the environment, i.e., rj/cu 2 . 3> l/O c - In fact, 
it turns out that inertia effects (finite capacitance) and 
a more refined modeling of the electromagnetic environ- 
ment lead for sufficiently large friction to only minor de- 
viations from the classical Smoluchowski predictio n 12 ' 13 
(they are relevant for a detailed quantitative analysis of 
actual circuits though) . For the quantum case considered 
in the sequel, the same is true if ft/3 > l/f2 c with H/3 be- 
ing at low temperatures the relevant scale for the coarse 
graining in timei^. 



III. CURRENT- VOLTAGE CHARACTERISTICS 

Mean values of relevant observables arc determined by 
the distribution P(6,t) determined from the QSE. Since 
most of the results can only be obtained numerically, 
we switch in this section to dimensionless quantities and 
scale energies in units of A, frequencies in units of A /h, 
and times in units of h/(Ap). In particular, this means 
to measure temperature in units of A./ks and currents in 
units of 7 C = Ej/Qq. The dimensionless QSE then reads 
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which is in fact a continuity equation for the probability, 
i.e., dP/Ot + OJ/06 = with the probability flux 
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Now, these expressions determine mean values with 
respect to phase (...) and time < ... >. For the current 
one has 



(HO)) 
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dtI(0)P(0,t) 
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and the voltage across the contact (V) = (9) follows as 

p27Z poo 

W)= d9 dtJ(6,t). (14) 

JO J-oc 

Due to the periodicity of the potential U(8,i) in phase 
and time [cf. Eq. ((3J)], one expands density and current 
according to 



n.k&Z 



J(6,t)= J2 J " 



ikQ-\-inujt 



(15) 



(16) 



n.keZ 



4 



with the normalization condition 



so that ((20)) takes the compact form 



P n .o = Sruo/Zn. 



Further, one writes due to 



1(9) = J2 Im{6,T)sm{m6) 



as well as 



D Q {0) = Dr. 



(9, t) cos(md) 



(17) 



(18) 
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This way, the QSE (|lip is cast in an algebraic equation 
for the expansion coefficients, namely, 

—LdPnk = hPn,k — i-^-(Pn-l,k + Pn+l.k) 
k Zp 

oo 

^ ^ ^7n (-^"n,fe— ni ^~n,fc+m) 
m— 1 

oo 

+ik D rn ,{P n , k ) (20) 

m ; — 

Practically, one works on a two-dimensional grid for n, k 
with | n | < -/V max , |fc| < K max . As already pointed out in^ 
the corresponding set of N max x K max coupled equations 
can be associated with a non-Hcrmitian lattice model for 
particles on a square lattice. In particular, one observes 
that there is a coupling between chains n and nil pro- 
portional to V ac , a coupling between chains k and k ± m 
proportional to the m-th harmonic of the Josephson cur- 
rent, and a coupling between chains k and k ± m! pro- 
portional to the m'-th harmonic of the quantum diffusion 
coefficient. Note that this latter coupling is absent in the 
classical regime. 

Now, the orthogonality of circular functions allows to 
express the mean current and voltage as 



< 1(9) > = J2 P o,fc J -fc 
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meaning that we only need to calculate Po,fc explicitly. 
In fact, peaks in this probability coefficient are related 
to the observed Shapiro steps of order n/k in the I — V 
characteristics. 

To solve (|2"0"|) numerically we define vectors 

Pn = (■ ■ ■ Pn,k, ■ ■ ■ , Pn,l, Pn,-1, ■ ■ ■ , Pn,-k, ■ ■ •) an d 

I = (...Ik,---, I\,I-x, I-k, ■ ■ ■) an d matrices 
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n )k,k' 
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-j-uP n ,k + hj h,k' 

Im (5k\k-rn ~ <5fe' ,k+m) 

ikD m > (6k',k-m' + h',k+m') (23) 



Ln^n = + ^ n+l) + S n , ~t . (24) 

2p 

This equation is solved via a recursive procedure (con- 
tinued fraction method-upward iteration) by introducing 

for n > the auxiliary quantity S n +iP n = ^(P n +i) 
and for n < by defining with n = —n the quan- 
tity Sn +1 ~^n = ^(Pn+i)- Accordingly, one has a 
simple equation for the relevant probability coefficients 
= [Lq - S x - 5 T ] + ~t, where 
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with the abbreviation p = (Vac/^p) 2 



IV. RESULTS 

We start with a brief discussion of actual experimen- 
tal parameters and proceed with a presentation of the 
numerical results. 



A. Approximations and parameters ranges 

We take typical experimental values for atomic con- 
tacts with Al electrodes^ with a superconducting gap 
A ~ 180peV. Temperatures are varied between about 
lOmK and lOOmK and the circuit is assumed to have an 
ohmic resistance of 200O such that p <C 1 and a capaci- 
tance on the order of fF. Typical microwave frequencies 
are Huj ~ 10~ 2 A — 10° A with u> > pA/h to observe frac- 
tional Shapiro steps. Within this range of parameters 
the conditions in (fTCTj) are fulfilled and phase diffusion is 
supposed to be affected by quantum fluctuations in the 
strong damping regime. In particular, /3E c /irp > 1 so 
that the energy scale related to friction hrj by far exceeds 
the thermal energy scale k^T. 

We now solve numerically the recursion (|25[) where 
convergence depends on the maximum number of spa- 
tial and temporal harmonics, the temperature and the 
external voltage. The numerics is quite sensitive at low 
temperatures and high voltages, but in the chosen ranges 
of these parameters accurate data can be achieved with 
iVmax = 90 and K max = 45. 



B. Numerical Results 

To analyze the Shapiro step structure we ca lculate nu- 
merically, from (f2"Tj) and (|2"2"j) . mean currents (1(9)) and 
mean voltages (V) (in the figures I and V respectively), 
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FIG. 2: I — V curves for an ac-driven atomic point contact 
with r = 0.995 and T = 0.005. Other parameters are w = 
2tt.10 -3 , 7? = 5, Vac = 5.10~ 3 , R = 10~ 3 (dimensionless 
units, see beginning of Sec lIII[) , 



over the temperature range specified above and for var- 
ious transmission coefficients. Figure [2] shows a typical 
I — V curve with the first integer and several fractional 
resonances in the quantum regime (low temperatures) 
and for a high transmissive junction. While resonances 
appear in a pattern very similar to the known classi- 
cal ones, the role of quantum fluctuations is revealed 
when one compares low temperature results obtained 
with the SE and those gained with the QSE, respectively 
(Fig. [3]). At higher temperatures the diffusion coefficient 
Dq —> Dei such that both descriptions deliver identical 
data, but there are substantial deviations at low tempera- 
tures where Dq(8) > Dei- Indeed, the classical equation 
predicts much sharper resonances than the quantum one, 
where the reduction in height and the increase in width 
is more striking at higher transmissions. This smearing 
out is basically absent away from the resonances. 

In order to have a better insight in this behavior 
the heights of the fractional peak Ii/ 2 [calculated from 
(A/2, max — A/2.min)/2 as differences between maximal 
and minimal peak values] are plotted as functions of tem- 
perature in Fig. [4] As already discussed, quantum fluc- 
tuations reduce the peak heights at lower temperatures 
and for higher transmissive channels. Since the overall 
peak structure is not altered, one may mislcadingly de- 
scribe a reduced height within the classical approach by 
an effectively enhanced temperature. However, while it 
is true that the quantum diffusion coefficient is typically 
larger than the classical one (Dq > Dei = ksT), due 
to its dependence on the phase 9 the QSE can in gen- 
eral not simply be reduced to the SE by replacing T by 
an effective temperature (see next section). Variations of 
peak heights with increasing transmission arc illustrated 
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FIG. 3: Fractional Shapiro steps for a high (r = 0.995, top 
panel) and lower (r = 0.99, bottom panel) transmissive chan- 
nels. Red(green)-solid lines depict data obtained with the 
QSE at T = 0.006(0.02), while blue(black)-dotted lines de- 
scribe results from the SE at T = 0.006(0.02). At a higher 
temperature T = 0.02 both approaches give identical curves. 
Other parameters are as in Fig. [2] 



in Fig. [5] for several fractional resonances. Interestingly, 
mean values I n lk saturate in the quantum case towards 
the ballistic limit t — > 1 with larger deviations from the 
classical data for higher order steps. 

Apart from reduced heights quantum effects appear 
as a widening of the resonances. A natural magni- 
tude to quantify this, is the full width at half max- 
imun (FWHM) of the peaks. As expected, we see in 
Fig. [5] that the spreading of the quantum peaks exceeds 
that of the classical ones at lower temperatures with the 
FWHM taking larger values for higher harmonics. Both 
predictions coincide only at relatively elevated tempera- 
tures. The relative strength of quantum fluctuations is 
larger at lower fractional steps, cf. Fig. [6] 

The crucial question is, of course, whether the influence 
of quantum fluctuations seen above could actually be de- 
tected in a real experimental set-up. In fact, for contacts 
with different transmissions also experimental data (see 
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FIG. 4: Resonant peak height vs. temperature for the Shapiro 
step 1/2 according to the classical approach (dotted) and the 
quantum one (solid) and transmission coefficients r = 0.995 
(shifted upwards by 0.01), r = 0.99, and r = 0.9. Other 
parameters are the same than in previous figures. 
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FIG. 5: Peak heights I„/ k 



1-T for resonances with 



n/fc=l/2, 1/3 and 2/3 (from top to bottom) at temperature 
T = 0.006. Solid lines correspond to the quantum case and 
dotted lines to the classical one. Parameters are the same 
than in previous figures. 



e.g. Ref fi£) deviate from results of the adiabatic classical 
theory. There are at least three possible explanations for 
this discrepancy: Landau-Zener (LZ) transitions between 
adiabatic surfaces, charging effects and associated quan- 
tum fluctuations, and spurious noise. With respect to the 
first one, it was shown in Refill that nonadiabatic LZ 
transitions enhance the magnitude of the supcrcurrcnt 
peak in almost ballistic channels. This effect is stronger 
at slightly elevated temperatures and in highly transmis- 
sive channels. Physically, nonadiabatic transitions be- 
tween £L and E + surfaces only occur if the diffusive 
passage of the phase through the LZ-range around 6 = tt 
is sufficiently fast compared to the instantaneous relax- 
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FIG. 6: FWHM vs. T for r = 0.995 and various fractional 
resonances (top) together with the relative strength of quan- 
tum fluctuations (FWHMq - FWHMci) / FW HMci (bot- 
tom). Other parameters are the same as in previous figures. 



ation time of momentum. Accordingly, LZ-transitions 
are suppressed towards very low temperatures in con- 
trast to what is observed for quantum fluctuations (see 
previous figures). Hence, we are left with either quan- 
tum fluctuations or spurious noise or, what is most likely, 
both to explain the addressed differences. Spurious noise 
alone can be captured by an effective temperature which 
is indeed the strategy that has been followed ir>i2. For 
this purpose, we analyze in the following to what extent 
this concept may also be applicable to effectively include 
quantum noise. 



V. QUANTUM DIFFUSION AND EFFECTIVE 
TEMPERATURE 



To describe the dynamics of the phase in the vicin- 
ity of a resonance, we consider instead of the current 
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biased circuit in Fig. Q] (Norton representation) the com- 
pletely equivalent circuit where the contact is voltage bi- 
ased (Thevenin representation). Within the adiabatic 
approximation and in the overdamped quantum range 
(3E c /irp >• 1 the voltage across the contact is then given 
by (again in physical dimensions) 

V = $ 6(t) = RI(0(t)) +V b + KcCosM) + RZ Q (9(t)) 

(26) 

with the voltage bias V b = IbR and the quantum noise 
Z Q {0) = y/pD Q {0)Zci. _ 

For a perfect voltage bias (no noise) and in absence of 
external driving the phase evolves as 9(t) = 9(0) + u>ot 
with the Joscphson frequency luq — 2eVb/h. In pres- 
ence of an ac-drive the n/fc-Shapiro resonance appears 
if kuio = nu at a corresponding voltage Vb = (V) = 
(n/k)<S>QLiJ. Thus, right at the center of the n/k resonance 
no dc-current flows and according to (|2l)|) the diffusive 
motion of the phase can be expressed as 



6(t) = 2z/sin(wt) + -u)t + S(t) 
k 



(27) 



with S being the stochastic component of the phase 
on top of the dominating deterministic part 6o(t) = 
2^sin(wi) + jwt where v = eV ac /hLo. Upon inserting 
this expression in (|26[) one arrives at 

71 

$ S(t) = RI(6(t)) +V b - -<S> q uj + Z Q (9(tj) . (28) 

Apparently, the dynamics of the stochastic part 5 is 
much slower than that of the deterministic part since 
\Vb — (n/k)ui\ <C (n/k)u> and 1(0) ps 0. This separation of 
time scales can be exploited when calculating time aver- 
aged currents. Namely, plugging the result ([27]) into the 
Fourier expansion (fT8|) of the current leads first to 



It it (r)< sm[ma n k(t)] \Jo(2mi>) 



m— 1 



+2 J2 P (2mv) cos(2pujt) + 2 cos[ma n k(t)] 
P >i 

x E J 2 P +i (2mv) cos[(2p + l)ut] > (29) 
p>i J 

with the abbreviation a n k(t) = (n/k)ujt + S(t) and the 
Bessel function J p . Taking now time averages over one 
period of the external drive and accounting for the time 
scale separation between deterministic and stochastic dy- 
namics of the phase we obtain 

I(9(t)) ps 2 £ Im{r) sm[m5(t)} 

m>l 

x J2 P (2mv)(cos (^—j—ujtj cos(2puit)) 
p>i 

-2 £ Im(r) sm[m5(t)} 



x £ ^2 P (2TO^)(sin i^Y^t) sin(2po;i))(30) 

p>0 



Here, the averages can only take two values, namely, 
or 1/2 depending on the relation between ^,2^ and 
(2p+l). Eventually, this yields the current phase relation 
in the vicinity of the n/k Shapiro resonance 

I(9(t)) ps E(- 1 )'" / ^( T ) JlnWkv) sin[ty fc (t)] , (31) 

Z>1 

where we put m = I fc, 2p = In with I £ N and intro- 
duced the scaled phase 4>k(t) = kS(t). 

The expression (j3Tj) is then inserted into (|28|) to pro- 
vide an approximate equation of motion for (f>k ■ While in 
general solutions are accessible only numerically, insight 
is already gained by keeping just the term with I = 1 in 
(EED, i-e., 



$o<fe = k [Vb — -$ w ) + R (-\) n kIkJn(2kv) suu> fc ) 



+kRZ Q (0 o + k /k) 



(32) 



In contrast to the classical case, here the noise term 
also depends on the phase. In the classical regime, 
one shows that (j32]) is identical to the equation of mo- 
tion for the phase in absence of ac-driving if parameters 
are renormalizecU^: I c ,&s = \klk J n (2kv)\, T ft = k 2 T, 
Vb :C s = k(Vb — ^$ w). This way, one gains repli- 
cas of the Ivanchenko-Zil'bcrman expression for the dc- 
supercurrent around each n/k peak, namely, 



\lc,eg Jeff 



(33) 



with fiz(x,y) = lTn{I 1 - lxy (y)/I- iX y(y)} the modified 
Bessel function of first kind. Now, in the quantum 
regime in leading order we may put Zq(9) ps Zq(9q) = 
\J (3Dq(9o) Zqi such that a similar renormalization ap- 
plies, however, with a modified temperature scaling, i.e., 



! 9,eff 



k 2 (pD Q (9 Q (t)))T >T cti . 



(34) 



This effective temperature depends also on the dissipa- 
tion strength, the driving frequency and amplitude, and 
is a nonlinear function of the actual environmental tem- 
perature. We note that the actual experimental data^ 
do not follow the scaling of T e ff with k, but rather can 
only be described by a much higher effective tempera- 
ture T e ff,exp which even affects the integer Shapiro reso- 
nances and has been attributed to spurious noise in the 
circuitry. The above enhancement due to quantum fluc- 
tuations may partially contribute to T c ff oxp , however, is 
not able to completely account for the discrepancy be- 
tween Tcffoxp and T cS . 

Beyond the case for I = 1 progress is achieved by 
assuming \Af3U"(9)\ < 1 so that one may expand 
[cf. Eq. 0] PD Q (9) ps 1 + k[3I'(9) with 

I' (9) w5^i«(-l) I »/ IJk (T)Ji n (2i*i/)cos^ fc ). (35) 
l>i 

Thus, if contributions with sufficiently large I are rele- 
vant, one may no longer replace cos(l<f>k) ~ > 1 meaning 
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that quantum noise cannot be captured by a global effec- 
tive temperature. Instead, its phase dependence leads to 
a local "temperature" and n/k resonances are not sim- 
ply replicas of the supercurrent peak. To extract signa- 
tures of this breakdown of the universal scaling behav- 
ior (|34|). the residual spurious noise dominating T e ff,exp 
must be substantially reduced. We note that the findings 
of Grabert et. ali£ who studied the supercurrent phase 
diffusion in absence of driving in the tunnel limit and 
at low temperatures within a time-dependent perturba- 
tion theory, obtained also an extension of the Ivanchenko 
Zil'berman expression similar to (|33p but with an effec- 
tive Josephson energy. This result was later reproduced 
within the QSE-formulation^. 

VI. SUMMARY 

In summary, the impact of quantum fluctuations on 
fractional Shapiro resonances, a hallmark of a non- 
sinusoidal current phase relations, is analyzed for atomic 
point contacts with highly transmitting channels in the 
presence of microwave. Known experimental I — V 
rcsult a 12 i 15 exhibit substantial deviations when compared 
with predictions from a classical adiabatic theory. While 
one explanation has been the appearance of spurious 
noise in the circuit, here, we find that quantum fluctua- 
tions may give rise to a similar effect. Departures from 
the classical approach become relevant for highly trans- 
missive channels and for sufficiently low temperatures. 



An effective description of quantum noise in terms of an 
effective temperature only applies to contacts with al- 
most sinusodial current-phase relations which may offer 
a way to distinguish between classical noise and quantum 
fluctuations in high transmissive contacts. As a prereq- 
uisite, however, unspecific spurious noise sources in the 
circuitry must be under control so that contacts are em- 
bedded in heat baths with temperatures of about T ~ 40 
mK or below. 

The results that we have obtained correspond to Al 
point contacts with a superconducting gap of Aai ~ 
200/16]/. For this metal our model predicts that quantum 
fluctuations play a pronounced role for very low temper- 
atures (T < 40 mK). Landau-Zener transitions are neg- 
ligible if T <C Tlz ~ 0.5AaiA - b ~ 1 K (for transmissions 
r > 1 — 10~ 4 )2i Experimentally, signatures of quantum 
fluctuations are expected to be even more dominant for 
materials with larger superconducting gaps such as e.g. 
Nb with ANb ~ 3600/16^/ leading to a clear separation 
between Tlz ~ 9 K and the temperature range where 
typical experiments arc performed (between 30 mK and 
150 mK). 
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